LETTER TO THE EDITOR 



Entanglement and tensor product decomposition for 
two fermions 

P Caban, K Podlaski, J Rembieliriski, K A Smoliriski and 
Z Walczak 

Department of Theoretical Physics, University of Lodz, Pomorska 149/153, 
90-236 Lodz, Poland 

Abstract. The problem of the choice of tensor product decomposition in a 
system of two fermions with the help of Bogoliubov transformations of creation 
and annihilation operators is discussed. The set of physical states of the composite 
system is restricted by the superselection rule forbidding the superposition of 
fermions and bosons. It is shown that the Wootters concurrence is not the proper 
entanglement measure in this case. The explicit formula for the entanglement of 
formation is found. This formula shows that the entanglement of a given state 
depends on the tensor product decomposition of a Hilbert space. It is shown that 
the set of separable states is narrower than in the two-qubit case. Moreover, there 
exist states which are separable with respect to all tensor product decompositions 
of the Hilbert space. 
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Entanglement is the key notion of quantum information theory and plays a 
significant role in most of its applications. The entanglement of a physical system 
is always relative to a particular set of experimental capabilities (see, e.g. P 
which is connected with decompositions of the system into subsystems. From the 
theoretical point of view this is closely related to possible choices of the tensor product 
decomposition (TPD) of the Hilbert space of the system. As a consequence the 
following question arises: How much entangled is a given state with respect to a 
particular TPD? 

In the present paper we discuss the problem of the choices of TPD in a system of 
two fermions, neglecting their spatial degrees of freedom and modifying tensor product 
in the rings of operators because of anticommuting canonical variables. We show that 
TPDs are connected with each other by Bogoliubov transformations of creation and 
annihilation operators. We also study the behavior of the entanglement of the system 
under these transformations. An importance of such investigation can be illustrated 
for example by the fact that the Bogoliubov transformations used in derivation of 
the Unruh effect also lead to the change of entanglement jHJ. Different approach to 
the entanglement in the system of two identical fermions, based on the asymmetric 
decomposition of the algebra generated by Oj, a\ (i = 1, 2) into tensor product of two 
subalgebras was taken up in 0]. Some aspects of the entanglement for two-fermion 
system were also discussed in 
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The theory of entanglement can be seen as the general theory of states 
transformations that can be performed on multipartite systems, with restriction that 
only local operations and classical communications (LOCC) can be implemented UJ. 
For the same reason, it was expected that additional restrictions should lead to new 
interesting physical effects and applications. Recently, it has been shown that such a 
restriction can be given by a superselection rule (SSR) 0IE1- 

In this work we restrict the set of physical states of the composite system by the 
requirement that we prohibit superpositions of fermions and bosons. This leads us to 
the SSR that is a weaker restriction (i.e., it admits larger set of states) than SSR based 
on the conservation of the number of particles 0. Moreover, we find the entanglement 
of formation taking into account the restriction imposed by our SSR. 

Let us consider the Hilbert space TL = C 4 with an orthonormal basis 
{|m, n}} m , n= o,i. With this basis we associate the following two operators: 

ai = |0,0)(l,0|-|0,l)<l,l|, (la) 
a 2 = |0,0)(0,l| + |l,0)(l,l|. (16) 

One can easily check that these operators and their Hermitian conjugations fulfill the 
following relations: 

{a i ,a J } = 0, {",.",! <>,,. i,j = 1,2, (2) 

where {., .} stands for anticommutator. Operators a] generate all the basis vectors 
from the "vacuum state" |0,0) via the relations 

|l,0) = a||0,0), (3a) 

|0,1> = 4|0,0), (36) 

|1,1) = 44|0,0), (3c) 

while the vacuum is annihilated by a,, i.e. a* |0, 0) = 0, i = 1, 2. We use the occupation 
number basis, i.e., the Fock basis, not the so called "computational basis". Thus with 
every orthonormal basis we can associate some representation of the algebra On 
the other hand it is clear that Eqs. @ can be interpreted as canonical anticommutation 
relations for two-fermion system. 

Every two orthonormal bases in 7i are connected by some unitary transformation 
belonging to the group E/(4). In the ring of operators these changes of bases are 
related to Bogoliubov transformations of creation and annihilation operators which 
will be discussed later on. 

One can naively expect that, as in the bosonic case, the operators a\ and a 2 should 
have the form a ® id and id 0a, respectively, where a is an annihilation operator for 
single fermion acting in C 2 , that is 

a|0}=0, a |1> = [0> , (4a) 
a t|0)=|l), a t |l}=0, (46) 
{a, a t } = id, (4c) 

and id denotes the identity operator. However it is not the case because a id 
and id ®a necessarily commute so they cannot fulfill the canonical anticommutation 
relations @ • In order to construct a± , a 2 out of the single annihilation operator a and 
to provide natural tensor product interpretation of basis vectors as \m, n) = \m) <g) |n) 
we have to modify only the tensor product of the operators acting in C 2 . Hereafter 
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we will denote the new tensor multiplication by the usual symbol ®. Such a modified 
tensor product is defined by the graded (supersymmetric) multiplication rule 

(A ® b)(o ® B) = (-l) F(a)F(b) Aa ® 65, (5) 

where a, b are monomials in a, aJ, i.e. a, b £ {id, a, cr, aa', a^a}, A, £? are arbitrary 
operators acting in C 2 and the "fermion number" F(a) is equal to the number of 
creation operators minus the number of annihilation operators building the monomial 
a, i.e. F(id) = 0, F(ot) = -F(a) = 1, F{aa^) = F(a^a) = 0. Consequently the 
Hermitian conjugation in this tensor product is of the form 

(a® b) f = (-l) F(a)F(b V® (6) 

The tensor multiplication introduced above is a special case of a more general 
structure known in mathematical physics as the braided tensor product [5]. As we 
can see from the relation (J5J, the new braided tensor product for monomials even in 
a, behaves like the standard tensor product. 

Finally, the relationship between the tensor product of operators and the tensor 



product of vectors is given by 

(id®id)(|m) ® \n)) = \m) <g> \n) , (7a) 

(a ® id)(|m) <g> \n)) = {-l) n a \m) ® \n) , (76) 

(a f ® id)(|m) ® |n)) = (-l)V |m) ® |n) , (7c) 

(id(8a)(|m) ® \n)) = |m) ®a|n>, (7d) 

(id®at)(|m) ® |n)) = |m) ® a f |n) . (7e) 

Now the annihilation and creation operators acting in the space C 2 ® C 2 and satisfying 
J2J take the desired form 

a\ = a ® id, <22 = id®a, (8a) 

= a 1 ' (g) id, a^id^a 1 . (86) 



Notice, that in the above equations ® denotes the new tensor multiplication thus 
Eqs. H4af) - l|7e| imply that operators l|8a|l -l|86 | fulfill the canonical anticommutation 
relations J5J). In particular, the matrix elements of operators (|8a|) - (|86| and Hla|) - ((T6t 
are identical in the basis n)} m n= o,i- 

Similarly, like in the case of quantum theory of fcrmionic fields in the system under 
consideration observables are restricted to combinations of even products of creation 
and annihilation operators. In particular the local observables are combinations of 
id(g>id and N\ = a^a ® id or id® id and N2 = id <g>a'a. It is implied by the SSR 
related to the requirement that the operator (— 1) F , where F is the fermion number 
operator, should commute with all observables |l(Jj . It means that superpositions of 
bosons and fermions are forbidden. In the basis (|3af) - (|Sc| (— l) F = diag{l, —1, — 1, 1}. 
Alternatively, this SSR is a consequence of the requirement that the squared time 
reflection operator must commute with all observables (see e.g. JI]). Indeed, 
antiunitary time inversion operator is defined here as follows 

TaiT -1 = a 2 , Ta 2 T _1 = — ai, (9) 
T|0,0) = |0,0). (10) 
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Thus T 2 = ( — 1) F . Due to the SSR the density matrix has to commute with (— 1) F , 
so the general state of this system is represented by the following density matrix 



(11) 




where w i ,v l > 0, £) i=1 ( 



1 and |6i| 2 < WiVi,i = 1,2. Consequently, possible 



states of subsystems obtained from by partial traces are 



Pi = 



P2 



wi+v 2 

U> 2 + V\ 

Wi + W2 

vi + v 2 



(12a) 



(126) 



Note that the diagonal form of H12al) - (|126|) is in conformance with the SSR in spaces 
of subsystems. Moreover, the states (I12al) - (|126l) exhaust all possible states of the 
subsystems. Therefore our subsystems are independent in the sense of the definition 
of the algebraic independence of subsystems This independence is due to the 

SSR (compare 4 where it was shown that in general algebras of observables of two 
identical fermions are nonindependent). The natural question arises: What is the form 
of the separable states for this system? According to Werner's definition the state 
is separable if it can be written in the form p = YliPiPi ® P21 where p\ and p\ are 
admissible states of subsystems and YliPi = 1' Pi — ®- Therefore, taking into account 
that p\ and p\ are of the form (|12a|) - (|12&j) . the separable states have the surprisingly 
simple diagonal form 



(13) 




with J2i Aj = 1, Aj > 0. Consequently, nondiagonal density matrices are nonseparable. 
Thus in this case the standard method of calculating entanglement measures should 
be taken with care. Indeed, as an example let us consider the Werner state 
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7 € [-1/3,1] 



(14) 



which belongs to the admissible states (|1 1|) . The Wootters concurrence ^H] of this 
state is equal to zero for 7 £ [—1/3, 1/3], therefore for two qubits the Werner state is 
separable for such values of 7. On the other hand, in our case this state is separable 
only when 7 = 0. Thus the Wootters concurrence does not define entanglement 
measure in our case. 

Instead, let us calculate directly the entanglement of formation [5], i.e.: 



E(p) = min^ piS(p\), 



(15) 



where S(pa) = ~ Tr pa log 2 pa is the von Neumann entropy and the minimum 
is taken over all the possible realizations of the state p — YliPi bPi) (V^il with 
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p l A = Ttb (\il>i) Taking into account the special form of the density matrix 

Ullfl we can find the explicit formula for the entanglement of formation 



E(p) = y^X w i + V i)$i 



(16) 



where 



S, 



and 



if' 



Vi and b j = 



(i-&)io g2 ^ + (i+e t )iog 2 ^ 



otherwise 



(17) 



6 = 



(18) 



if 7 ^ 0, 



(19) 



yj{wi - Vi) 2 +4|6i, 

It is interesting that the formula similar to l|17|) was obtained in f° r the so called 
correlational entropy of the two-level system. Note that maximal value of E(p) is 
equal to 1. In the case of the Werner state 114|) the entanglement of formation (|16|) is 

r 1+7 

E(p w ) = I 2 

[ if 7 = 0. 

Thus, as expected, E(p\y) ^ for entangled (nondiagonal) states and E{pw) = 
for a separable (diagonal) state. For 7 = 1 we have the maximally entangled Werner 
state. Notice, that the restriction of admissible states by SSR implies that in our case 
we have no asymmetry in the definition of the entanglement of formation, in contrast 
to observations of 0]. 

Let us consider the problem of decomposition of our system into two subsystems. 
Such a decomposition corresponds to the different choices of canonical variables cii, a\. 
This is extremely important because each choice of et^ , a\ defines in the Hilbert space 
H the corresponding tensor product structure (|^|l- (f7ejl such that the creation and 
annihilation operators take the form analogous to (|8o|l (|8b|l . Each TPD defines a set 
of local observables of the form A <g> id and id(g)B [cf. the discussion after H8a |l - I|8fe|l ]. 
Moreover, the notion of a local observer is determined by his experimental access to 
local observables (see e.g. £Q). 

Different choices of canonical variables ai , a\ are connected by transformations 
which preserve the canonical anticommutation relations J5J (Bogoliubov transforma- 
tions!). Therefore Bogoliubov transformations give us all possible decompositions of 
the two-fermion system into two subsystems (two fermions). Such decompositions of 
the system correspond to the tensor product decompositions of the space Ji = <C 2 Cg)C 2 , 
appropriate to the definition of the subsystems. In the case under consideration the 
problem of finding all possible TPDs consistent with our SSR is equivalent to deter- 
mining all possible Bogoliubov transformations commuting with the superselection 
operator T 2 = (— 1) F . 

Let us notice first that operators ai,a\ in every orthonormal basis can be 
represented in the form (J3^— and vice versa such operators define orthonormal 



f By Bogoliubov transformations we mean here all transformations of creation and annihilation 
operators (linear as well as nonlinear) which do not change the canonical anticommutation relations. 
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basis via (|3a|) - (|3"c| . Thus different choices of these operators are connected with 
different choices of orthonormal bases in the Hilbert space. Therefore 

a- = U ai U\ (20) 

where U is an unitary matrix. As we have mentioned above, the consistency with 
SSR means that U commutes with T 2 = (— 1) F = diag{l, — 1, — 1, 1}. So U can be 
represented as the following product of unitary matrices 

/( -w* 



U = 





(21) 



10 
1 

\w o c* 

with \a\ 2 + \/3\ 2 — 1 and |£| 2 + M 2 = 1, where we took into account the fact that 
Eq. (|20|l determines U up to an overall phase. Thus all Bogoliubov transformations 
admissible by the SSR form the group SU(2) <S> U(2). Applying these transformations 
to the explicit matrix form of aj,aj calculated from (|lqf) - (|lb|) . one can show that in 
the ring of creation and annihilation operators the transformations (|20|l are realized 
as 



SU (2) transformations which do not mix creation and annihilation operators 

(22) 



a 'i \ _ ( a P \ ( fl i 
a'J \-0* a* J {a 2 

SU (2) transformations which mix creation and annihilation operators 

a[ \ _ f C uj \ f ct\ 



A ) V " w * C* J V °2 



(23) 



h 2 

• nonlinear one-parameter transformations 

a[ = a x e 1 ^ = oi[l + (e** - 1)N 2 ], (24a) 
a' 2 = a 2 e lxNl = a 2 [l + (e** - l)N{\. (246) 

Note that (|24a|) (|246|l for x = 7r are the so called Klein- Wigner transformations (cf. 
0]). The Bogoliubov transformations which lead to physically distinguishable TPDs 
should change the local observables 2Vj = ajflj. Therefore such transformations have 
the form: (22J with both a ^ 0, /3 ^ and/or (J^ with both C ^ 0, w ^ 0. 

Now, the natural question arises: How does the same state look like for local 
observers connected with different TPDs? The answer is quite obvious: If their TPDs 
are connected by Bogoliubov transformations then density matrices representing the 
state are connected by similarity transformations, i.e. p' = U pW . However, in general, 
such transformation change the entanglement measure E(p), i.e. entanglement depends 
on the choice of TPD (and hence the local observers). In particular, for any state, there 
exists a pair of observers for whom this state is separable, since the density matrix (fill) 
can always be diagonalizcd by means of the transformations 121f) . We point out that 
there exists a class of supers eparable states p ss = ^diag{s, 1 — s, 1 — s, s}, s G [0, 1], 
which are separable for every observers. Note also that in the case of two qubits only 
one superseparable state exists, namely the maximally mixed state po = ^1. 

Now we show that it is possible to construct dynamics consistent with our SSR. 
For such a dynamics admissible TPDs are related to symmetries of the Hamiltonian. 
An example of that dynamics is the Thirring model |17| in 1+0 dimensional space-time 
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describing a fermionic quantum mechanical system. The corresponding Lagrangian is 
of the form: 

2 2 

l = YtWlWi - m 4^) - a(E ^M- ( 25 ) 

i=l i=l 

The solutions of the equations of motion derived from the Lagrangian l|25|) are 

^(t)^a t e-^ m+x+2XN ^\ i+j, (26) 

where the time-independent operators a; and a\ satisfying J2J can be represented in 
the form H8a fl -Q8fr [l . The Hamiltonian of this system is 

H=(m + \)(Ni+N 2 ) + 2XN 1 N 2 (27) 

and describes two fermionic oscillators with the quartic interaction term. Notice that 
T 2 = (— 1) F commutes with H, thus Thirring model dynamics undergoes our SSR. 

In the special case of A = —\m all the Bogoliubov transformations i|22|) - (l24b(l 
form the symmetry group of H, i.e. H(a,a' i ) = H(a',a ,Jf ). Thus, this symmetry 
group gives us a freedom with a choice of a concrete decomposition of the system into 
two subsystems. Consequently the related TPDs are connected by the Bogoliubov 
transformations (j2"2"|l - lj24&l) . 

In conclusion, we have investigated the dependence of entanglement for two- 
fermion system on tensor product decompositions in the presence of the superselection 
rule. We have shown that the Wootters concurrence is not a proper entanglement 
measure in this case. The crucial point in finding an explicit form of entanglement of 
formation for such a system was determining the states of subsystems, admissible by 
the superselection rule. We would like to stress that these states are not qubit states. It 
is interesting that the set of separable states is narrower than in two-qubit case, namely 
it consists of only the states represented by diagonal density matrices. Moreover, we 
found the class of superseparable states, i.e. the states which are separable with respect 
to all tensor product decompositions of Hilbert space. 
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